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a-cluster correlations in the ground states of 12 C and 16 O are studied. Because of the a cor- 
relations, the intrinsic states of 12 C and 16 have triangle and tetrahedral shapes, respectively. 
The deformations are regarded as spontaneous symmetry breaking of rotational invariance, and the 
resultant oscillating surface density is associated with a density wave (DW) state caused by the 
instability of Fermi surface with respect to a kind of lp-lh correlations. To discuss the symmetry 
■ breaking between uniform density states and the oscillating density state, a schematic model of a 

few clusters on a Fermi gas core in a one-dimensional finite box was introduced. The model analysis 
suggests structure transitions from a Fermi gas state to a DW-like state via a BCS-like state, and to 
(N. a Bose Einstein condensation (BEC)-like state depending on the cluster size relative to the box size. 

{^jqj It was found that the oscillating density in the DW-like state originates in Pauli blocking effects. 
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I. INTRODUCTION 



Nuclear deformation is one of the typical collective motions in nuclear systems. It is known that ground states of 
nuclei often have static deformations in the intrinsic states, which are regarded as spontaneous symmetry breaking of 
the rotational invariance due to collective correlations. Needless to say, the broken symmetry in the intrinsic states is 
restored in nuclear energy levels, because total angular momenta are good quanta in energy eigenstates of the finite 
£jj 1 system. Not only normal deformations of axial symmetric quadrupole deformations but also triaxial, octupole, and 
super deformations have been attracting interests in these decades. To investigate those deformation phenomena 
C ■ mean-field approaches have been applied, in particular, for heavy mass nuclei. 

In light nuclear systems, further exotic shapes due to cluster structures have been suggested. For instance, a triangle 
shape in 12 C and a tetrahedral one in 16 have been discussed based on the cluster picture that 12 C and 16 are 
£> | considered to be 3a and 4a systems. In old days, to understand spectra of 12 C and 16 non-microscopic a-cluster 
models have been applied P, H[ . From vibrations of the triangle structure of thee a particles and the tetrahedral one of 
four as, Wheeler has suggested the low-lying J = 3 states in 12 C and 16 pj, which are now considered to correspond 
to the lowest negative-parity states 12 C(3j~, 9.64 MeV) and 16 0(3 1 ', 6.13 MeV) established experimentally. In 1970's, 
cluster structures of the ground and excited states in 12 C and 16 have been investigated by using microscopic and 
semi- microscopic^ cluster models , a molecular orbital model [l2| , and also a hybrid model of shell model and 



cluster model [13j. 

For 12 C, the ground state is considered to have the triangle deformation because of the 3a-cluster structure. In 
addition, a further prominent triangle 3a structure has been suggested in 12 C(3j~, 9.64 MeV). Although the cluster 
structure of the ground state, 12 C(0^), may not be so prominent as that of the 12 C(3 1 "), the J w = 0^ and 3^ states 
are often described by the rotation of the equilateral triangle 3a configuration having the symmetry. In contrast 
to the geometric configuration suggested in 12 0(0]*") and 12 C(3j~), a developed 3a-cluster structure with no geometric 
configuration has been suggested in the 0^ state assigned to 12 C(0^, 7.66 MeV) by (semi-)microscopic three-body 
calculations of a clusters H Gil- I n the 0^" state, three a particles are weakly interacting like a gas, for which 
the normal concept of nuclear deformation may be no longer valid. For the 0^ state, Tohsaki et al. has proposed a 
new interpretation of a dilute cluster gas state, where a particles behave as bosonic particles and condensate in the 
S orbit [14[. This state is now attracting a great interest in relation with the Bose Einstein condensation (BEC) in 
nuclear matter [l5j . 

Let us consider the cluster phenomena in 12 C from the viewpoint of symmetry breaking. If the symmetry of the 
rotational invariance is not broken, a nucleus has a spherical shape. In the intrinsic state of 12 C(0f ), the spherical 
shape changes to the triangle shape via the oblate shape because of the a-cluster correlation. It is the symmetry 
breaking from the rotational symmetry to the axial symmetry, and to the Z?3h symmetry. In the group theory, it 
corresponds to 0(3) — > 0(2) x — > where the symmetry breaking from the continuous (rotational) group to 
the discrete (point) group occurs. In the excited state, 12 C(0^~), the system again may have the continuous group 
0(3) symmetry. It indicates that the cluster correlations in 12 C(0^) and 12 C(0^) may have different features in terms 
of symmetry breaking. The triangle shape with the symmetry in 12 C(0^) is characterized by the geometric 
configuration, while the 12 C(0^) has no geometric configuration. Now a question arises: what is the mechanism of the 
symmetry breaking of the continuous group in 12 C(0|" ), which is restored again in 12 C(0^). One of the key problems 



2 



is the geometric configuration because of a correlations in the ground state. The triangle state has oscillating surface 
density along the oblate edge. It can be understood by the density wave (DW)-like correlation caused by the lp-lh 
correlation carrying a finite momentum in analogy to the DW in infinite matter with inhomogeneous periodic densities, 
which has been an attractive subject in various field such as nuclear and hadron physics jl6l434| as well as condensed 
matter physics [35l. |36|. Indeed, in our previous work, we interpreted the triangle shape as the edge density wave on 
the oblate state by extending the DW concept to surface density oscillation of finite systems (3tJ • Then the structures 
of 12 C(0+) and 12 C(0^) may be associated with the DW and the BEC phases in infinite matter, respectively. The 
mechanism of the geometric triangle shape in the finite system may give a clue to understand an origin of DW in 
infinite matter. 

Similar to 12 C, a geometric configuration with a tetrahedral shape in 16 has been discussed in theoretical studies 
with cluster model calculations 0, El El and also with Hartree-Fock calculations [38l - l40j . The excited state, 16 0(3^, 
6.13 MeV), is understood by the tetrahedron vibration or the rotation of the tetrahedral deformation with the Td 
symmetry, while the static tetrahedral shape in the ground state has not been confirmed yet. The tetrahedron 
shape in 16 O(0+) and 16 0(3 1 ") is supported in analysis of experimental data such as E3 transition strengths for 



3j — » 0^ }4lJ_and a-transfcr cross sections on 12 C 42|. The tetrahedral shape tends to be favored in cluster-model 



calculations [3j, [11[. However, Hartree-Fock calculations with tetrahedral deformed mean-field potentials usually 
suggest the spherical p-shell closed state as the lowest solution for 16 O except for the calculations using effective 
interactions with specially strong exchange forces (38l - l40j . If the ground state of 16 O has the tetrahedral shape, it 
may suggest the breaking of the 0(3) symmetry into the Td symmetry. In the excited + states of 16 0, 12 C+a cluster 
structure was suggested in the 0^" state at 6.05 MeV 0, H, US S; 01 • Moreover, in analogy to the 3a-cluster gas state 
of 12 C(0^~), a 4a-cluster gas state with the a condensation feature has been suggested recently in a + state above 
the 4a threshold [HI, H^. Similarly to 12 C, the possible tetrahedral shape in 16 may lead to symmetry breaking of 
the continuous group in 16 O(0^), which is restored in higher + states. Again, the geometric configuration due to a 
correlations in the ground state should be one of the key problems. 

Our aim is to clarify the a-cluster correlations with geometric configurations in the ground states of 12 C and 
16 O, and understand the mechanism of the symmetry breaking from continuous (rotational) groups into discrete 
(point) groups. We first confirm the problem whether the triangle and tetrahedron shapes are favored in the intrinsic 
states of 12 C and 16 0. For this aim, we use a method of antisymmetrized molecular dynamics (AMD) j47l - l49j and 
perform microscopic many-body calculations without assuming existence of any clusters nor geometric configurations. 
Variation after the spin-parity projections (VAP) is performed in the AMD framework [5(3]. The AMD+VAP method 
has been proved to be useful to describe structures of light nuclei. With this method, one of the authors, Y. K-E., has 
succeeded to reproduce various properties of ground and excited states of 12 C [13, HU, and confirmed the formation 
and development of three a clusters in 12 C in the microscopic calculations with no cluster assumptions for the first 
time. The result was supported also by the work using the method of Fermionic molecular dynamics (H^j, which shows 
a similar method to the AMD. 

In this paper, we apply the AMD+VAP method to 16 as well as 12 C and analyze the intrinsic shapes of the 
ground states. We show that the geometric configurations having the approximate and Td symmetry arise in 
the ground states of 12 C and 16 0, respectively. To discuss appearance of the geometric configurations, we perform 
analysis using a simple cluster wave functions of Brink-Bloch (BB) a-cluster model [53|, while focusing on Pauli 
blocking effect on rotational motion of an a cluster. Important roles of the Pauli blocking effect in appearance of 
geometric configurations are described. We also introduce a schematic model by considering clusters on a Fermi gas 
core in a one-dimensional (ID) finite box, which can be linked with clusters at surface in a 3a system. By analyzing 
the ID-cluster wave function, in particular, looking at Pauli blocking effects from the core and those between clusters, 
we try to conjecture what conditions favor BCS-like, DW-likc, and BEC-likc correlations. 

This paper is organized as follows. In the next section, intrinsic shapes and cluster formation in the ground states 
of 12 C and 16 are investigated based on the AMD+VAP calculation. In Sec. Mil a Pauli blocking effect in a-clustcr 
systems and its role in a-cluster correlations is described by analysis of BB a-cluster wave functions. In Sec. IIV1 the 
schematic model of clusters on the Fermi gas core in the ID finite box is introduced and the roles of Pauli blocking 
effects in a-cluster correlations are discussed. Summary and outlook are given in Sec. |V] The relations between 
3a- and 4a-cluster wave functions and triangle and tetrahedral deformed mean-field wave functions are explained in 
appendix [X] In appendix IB"! and [Cl features of weak-coupling wave functions in the ID-cluster model are described. 



II. SHAPES AND CORRELATIONS IN THE GROUND STATES OF 12 C AND ls O 

We discuss here intrinsic deformations of the ground states of 12 C and 16 O based on the AMD+VAP calculation. 
The AMD method has been applied for various light mass nuclei and has been successful in describing cluster structures 
as well as shell- model-like structures in light- mass nuclei. In the present work, the AMD+VAP method, i.e., variation 
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after spin and parity projections in the AMD framework, is applied to 12 C and 16 0. For the details of the framework, 
the reader is refereed to, for instance, Refs. [49l. |50|. 

A. Variation after projection with AMD wave function 

In the AMD framework, we set a model space of wave functions and perform the energy variation to obtain the 
optimum solution in the model space. 

An AMD wave function is given by a Slater determinant of Gaussian wave packets, 

$amd(Z) = -=A{<pi,<p2, ¥>a}, (1) 
y/A\ 

where the ith single-particle wave function is written by a product of spatial ((/)), intrinsic spin and isospin wave 
functions (x T ) as 

w = faxlx!, (2) 

2„\ 4/3 , , X,- 



M r i)=[-) cxpj-Kr,-^) 2 }, (3) 

xl = (\ + Si)x f + {\-iU)xi r - ( 4 ) 

^>X; and xf are spatial and spin functions, and \\ is the isospin function fixed to be up (proton) or down (neutron). Ac- 
cordingly, an AMD wave function is expressed by a set of variational parameters, Z = {Xi, X2, • • • , Xyt, £1, £2, • • • , 
The width parameter v relates to the size parameter b as v = 1/26 2 and it is chosen to be v = 0.19 fm that minimizes 
energies of 12 C and 16 0. 

The center positions Xi, X2, • • • , X^ of single- nucleon wave packets are independently treated as variational param- 
eters. Thus existence of any clusters are not a priori assumed in the AMD framework. Despite of it, the model wave 
function can describe shell-model structures and cluster structures because of the antisymmctrizer and the flexibility 
of spatial configurations of Gaussian centers. If a cluster structure is favored in a system, the corresponding cluster 
structure is automatically obtained in the energy variation. 

For even-even nuclei, the ground states are known to be J n = + states, i.e., they are symmetric for rotation. 
Intrinsic deformation is understood as spontaneous symmetry breaking with respect to the rotational invariance, 
which is restored in the J T = + ground states. It means that when an intrinsic state has a deformation the ground 
state is constructed by the spin and parity projections from the intrinsic state. In more general, spin and parity 
are good quanta in energy eigenstates of nuclei because of the invariance of the Hamiltonian for rotation and parity 
transformation. Therefore, to express a J T state, an AMD wave function is projected onto the spin-parity eigenstate, 

<J>(Z) = P// a -$amd(Z), (5) 

where P{j K is the spin-parity projection operator. 

To obtain the wave function for a J™ state, the energy variation is performed for the spin-parity projected AMD wave 
function $(Z) with respect to variational parameters {Z}. This method is called variation after projection (VAP). The 
AMD+VAP method has been applied to various light nuclei for structure study of ground and excited states. For the 
ground states of 12 C and 16 0, we perform the variation of the energy expectation value ($(Z)|i?|$(Z))/($(Z)|$(Z)) 
for the J™ = + projected wave function and get the optimum parameter set {Z} that minimizes the energy. Then, 
the AMD wave function $amd(Z) given by the optimized {Z} is regarded as the intrinsic wave functions of the ground 
state. 

An AMD wave function is expressed by a single Slater determinant; however, the spin-parity projected wave function 
is no longer a Slater determinant but it is a linear combination of Slater determinants except for the case that the AMD 
wave function before the projection is already a spin-parity J n eigenstate. If the intrinsic state has a deformation the 
projected wave function contains some kinds of correlations beyond Hartree-Fock approximation. 

B. Intrinsic structures of 12 C and 16 O 

We apply the AMD+VAP method to the ground states of 12 C and 16 O, and discuss their intrinsic structures. 
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1. Density distribution 



The ground states of 12 C and 16 have the intrinsic deformations. The density distribution of the intrinsic wave 
functions 3>amd are shown in Fig. [T] The result for 12 C shows a triaxial deformation with a triangle feature, while 
16 has a deformation with a tctrahcdral feature. The quadrupolc deformation parameters (/3,j) evaluated by the 
quadrupole moments are (/3, 7) = (0.31,0.13) for 12 C and (/?, 7) = (0.25,0.09) for 16 0. The triangle and tetrahedral 
shapes are caused by a-cluster correlations. Strictly speaking a clusters are not ideal (0s) 4 clusters but somewhat 
dissociated ones. Moreover, the triangle and tetrahedron arc not regular but distorted as an a cluster is situated 
slightly far from the other as. 




FIG. 1: (color on-line) Density distributions in the intrinsic states for the ground states of (top) 12 C and (bottom) ie O obtained 
by the AMD+VAP calculation. The density integrated on the z, x, and y axes is plotted on the x-y, y-z, and z-x planes. 



2. Oscillation of the surface density 

Let us consider the deformations of 12 C and 16 from the viewpoint of symmetry breaking. The highest symmetry 
is the sphere, which is realized in the closed p 3 / 2 -shell and p-shell states for 12 C and 16 0, respectively. As a higher 
symmetry can break into lower symmetry due to multi-nucleon correlations, the deformation mechanism of 12 C is 
interpreted as follows. Because of the a-cluster correlations, the rotational symmetry in the sphere breaks to the 
axial symmetry in an oblate deformation and changes into the D^h symmetry in the regular triangle 3a configuration, 
and then it breaks to the distorted triangle. The symmetry change, spherical— >oblatc— >triangle, corresponds to 
0(3) —¥ 0(2) x Z2 — s- £>3h- Similarly, the deformation of 16 is understood as the rotational symmetry in the sphere 
breaks into the Td symmetry in the regular tetrahedral 4a configuration, and it breaks to the distorted tetrahedron. 
The structure change from the spherical to the tetrahedron is the breaking of the 0(3) symmetry to the Td symmetry. 
Note that the continuous (rotational) groups break to the discrete (point) groups in the triangle and tetrahedron 
deformations. 

We discuss the connection of cluster correlations with the density wave, which is characterized by the static density 
oscillation at the surface. The surface density oscillation occurs at the symmetry breaking from the axial symmetry 
(oblate shape) to the Z?3h symmetry (triangle) in 12 C and that from the rotational symmetry (sphere) to the Td 
symmetry (tetrahedron) in 16 0. The triangle deformation contains the (Kj — Y^ 3 )/\/2 component with the Y® 
deformation in the density, while the tetrahedral one has the (\/5Y® + s/2Y^~ — V2Y^~ 3 )/3 component, which can be 
transformed to (Kj -2 + Y 3 +2 )/v / 2 by the rotation. Indeed, as described in appendix [A] an ideal 3a(4a)-cluster wave 
function with the triangle (tetrahedral) configuration has the density having the finite components of (Y^~ 3 — Y^~ 3 ) / y/2 
((V0T3 + \/2Y^~ 3 — v / 2^3 +3 )/3) and it can be described by the DW-type particle-hole correlations in case of weak 
deformations. Note that the DW in the triangle shape is characterized by the particle- hole correlations on the Fermi 
surface carrying the finite angular momentum (I, m) = (3, ±3) and that in the tetrahedral one is given by the particle- 
hole correlations with (l,m) = (3, ±2) as described in the appendix. Namely, the symmetry breaking is characterized 
by the (Y^ 3 — Y^ 3 )/\/2 component whose amplitudes linearly relate to the order parameter of the DW correlations 
as shown in Eqs. (|A3|) . (| A6|) . and (|A7[) . In other words, the symmetry broken states have finite (Y3 -3 — Y^~ 3 )/\/2 
components in surface density and show the oscillation density with the wave number three. 
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To analyze the surface density oscillation in the ground states of 12 C and 16 O, we perform the multipolc decompo- 
sition of the intrinsic density obtained with the AMD+VAP calculation 

P (r = Ro, 9, 0) = p(R ) Y, a l™Yr(6, </>), (6) 
i in 

at a certain radius r = R 0: and discuss the (Y^~ 3 —Y^~ 3 ) / y/2 components. We take Ro to be the root mean square radius 
of the intrinsic state. p(i?o) is determined by the normalization ctoo — P and a im has the relation ai m = ( — l) m o:*_ m 
because the density p(r = Ro,9, 4>) is real. 

The density plot at r = Rq on the 9-(f> plane for 12 C and that at r = Ro and 9 = ir/2 as a function of <fi are shown in 
Fig. [3J As seen clearly, the density on the oblate edge shows the oscillation with the approximate wave number three 
periodicity, which comes from the a-cluster correlation. In the right panel of Fig. [5J we also plot the density for the 
ideal symmetry given only by the (I, m) = (0, 0) ,(2,0), (3, 3) and (3,-3) components to show somewhat distortion 
of the AMD+VAP result from the ideal symmetry. In the coefficients \ai m \ shown in Fig. @J it is found that the 
y ; ±3 components are actually finite indicating that the axial symmetry is broken to the triangle shape in 12 C. 

For the density in 16 0, the 9-(p plot is shown in Fig. [31 and the coefficients of the multipole decomposition are 
shown in Fig. [H The tetrahedral component y/5Y®/3 + v / 2Y 3 +3 /3 — v / 2Y 3 _3 /3 is shown by the hatched boxes at 0:30 
and «33 in Fig. [4] The open boxes indicate the distortion components from the tetrahedron. The distortion exists 
in the axial symmetry components, a%o, ct2o, and aio, coming from the spatial development of an a cluster from the 
others as explained before. 

Thus, the intrinsic states of 12 C and 16 show the surface density oscillation with the (Y^~ 3 — Y3 )/ V% component. 
The wave number three oscillation characterized by the (Y^ 3 — Y^ 3 )/\/2 component is understood by the a-cluster 
correlations with triangle and tetrahedral deformations, which are interpreted as the DWs on the oblate and spherical 
shapes, i.e., the spontaneous symmetry breaking of axial symmetry — > and rotational symmetry — > Td symmetry. 




FIG. 2: (color on-line) Left: density at r = Ro for 12 C calculated by the AMD+VAP is plotted on the 0-cj> plane. Right: that at 
r — Ro and 9 = 7r/2 line (the solid line). The density for the ideal Z?3h symmetry given only by the (l,m) = (0,0) ,(2,0), (3, 3) 
and (3,-3) components is also plotted (the dashed line). Ro is taken to be 2.53 fm. 




e (Ttrad) 



FIG. 3: (color on-line) density at r = Ro for 16 C calculated by the AMD+VAP is plotted on the 6-(f> plane. Tio is chosen to 
be 2.81 fm. 
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FIG. 4: The coefficients aim of the multipole decomposition of the intrinsic density of (top) 12 C and (bottom) 16 O calculated by 
the AMD+VAP. The hatched areas in the bottom panel are the tetrahedron component -JEY® /3+V2Y^~ 3 /3-\/2Yf 3 /3 defined 
by O3o tch = \/572a33 and a- 1 !' 1 * = a 33 . The open area for the Y§ component is defined by the relation a 30 = a^ cn + a3o tch . 



III. PAULI BLOCKING EFFECT IN a CORRELATIONS 



As shown in the previous section, the intrinsic states of 12 C and 16 obtained with the AMD+VAP calculation 
contain the a-cluster correlations with the triangle and tetrahedral deformations. For 12 C, the triangle deformation 
has been suggested in 3a-cluster models, in which 3a state with a regular triangle configuration is the energy minimum 
state [1,0, EH- I n contrast to such the geometric configuration as the triangle, the second + state of 12 C is considered 
to be a cluster gas state where 3 a clusters are freely moving in dilute density like a gas without any geometric 
correlations between clusters @, H, H3, EII • The J state is associated with the a condensation in analogy to BEC 



because it is regarded as a system of three as occupying the same S orbit [14 1 . 

From the viewpoint of symmetry breaking, the symmetry is broken in the state with the triangle deformation 
and it is restored in the 0^" state. In the 0^ state, the axial symmetry breaks down to the D^h symmetry due to the 
3a-cluster structure. One of the characteristics of the 0+ state is the oscillating surface density caused by the angular 
correlation of a clusters. Then the structure change from the Of to the 0^" is expected to connect with the transition 
from the symmetry broken state with the angular correlation to the symmetric state with no (or less) correlation 
between a clusters. 

The origin of the angular correlation in the 0^~ state and the transition into the uncorrelated 0^" state can be 
understood by the Pauli blocking effect between clusters as follows. Let us here consider motion of an a cluster 
around a 2a core in the BB 3a-cluster model wave function. In the BB 3a-cluster model, a clusters are located 
around certain positions Si, S2, and S3, and the wave function is written as 

$bb(Si, S 2 , S 3 ) = —=A{U ra (t) Sl X T< j(j)s. } X Ta (t>s 3 X T(7 } , (7) 

where X Ta is the spin-isospin wave function with r = {p,n} and a = {1\4-}- We assume that 2 a clusters placed at 

51 = (0, 0, d/2) and S 2 = (0, 0, —d/2) form the core. The third a is placed at S3 = (0, y, z) for y = rcos9 y , z = rsm9 y 
on the (y, z)-planc (sec Fig. [5]). Because of Pauli blocking effect between clusters, the motion of the third a around 
the core is restricted in the Pauli allowed region. Particularly when the a exists near the core, rotational motion is 
strongly blocked because of the existence of other a clusters. The Pauli allowed and forbidden areas for the rotation 
of the angle 9 y for the third a center in the (y, ,z)-plane are presented in Fig. [5] In the figure, the norm Nsaiu, z ) = 
($ BB (Si, S 2 , S 3 )|$sb(Si, S 2 , S 3 )) of the BB wave function $bb(Si, S 2 , S3) with the parameters, Si = (0,0, d/2), 

5 2 = (0, 0, —d/2), and S3 = (0, y, z) is shown in the {y, z)-plane. d is fixed and taken to be a small value. The norm 
is normalized by the value for the a on the y-axis, 

- Af 3a {y = r cos 6 y ,z = r sin 9 V ) 

A/ 3Q (y , z = --r— y - y — . (8 

N3a(y = r,z = 0) 

The area with J\f 3a ~ 1 is the allowed region where the a feels no Pauli blocking with respect to the rotational 
motion, while the ~ region is the blocked region where it feels the strong Pauli blocking effect from a clusters 
of the core. It means that, when the third a exists near the core, its angular motion is blocked by the 2a core. As 
a result, the third a is confined in the Pauli allowed region around the y axis, and it has the angular correlation 
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against the 2a direction. Consequently, a compact 3a state has a geometric structure of the triangle deformation 
and it has the surface density oscillation. On the other hand, as the cluster develops specially the Pauli blocking 
effect becomes weak. When the a is far enough from the core, it can freely move in the rotational mode. Then the 
angular correlation vanishes and the system transits to angular-uncorrelated state. We note that in cluster physics 
this transition is known to be the change between strong cluster coupling and weak cluster coupling states, where the 
angular momentum of the inter-cluster motion couples with inner spins of clusters strongly and weakly, respectively 
(see Fig. [5]). What we call the cluster coupling is coupling between clusters and it is different from the terminology of 
strong and weak couplings in the BEC-BCS crossover phenomena, which relate to the coupling between nucleons in 
a pair (or in a cluster). 

The Of state of 12 C is considered to be the compact 3a state in the strong Pauli blocking regime and corresponds 
to the angular correlated state with the surface density oscillation due to the triangle deformation. In the 12 C(0.f ), 
clusters spatially develop well and the system goes to the non-angular-correlation state. Strictly speaking, in the 
12 C(0^), all clusters develop to form a uncorrelating three a state associated with the a condensation, where the 
radial motion is also important as well as the angular motion. Nevertheless, we can say that, for the restoration of the 
broken symmetry with the surface density oscillation in the 0^ to the rotational symmetry in the 0^", the transition 
from the angular correlated state to the uncorrelated state in the cluster development is essential. 

In the above discussion, we consider the angular motion in the intrinsic (body-fixed) frame, i.e., the y-z plane. 
Since the system is symmetric for the rotation around the z-axis, the motion of the third a is free for the z-rotation. 
This rotational mode around the z axis is nothing but the projection onto the J z = 0, which is usually performed in 
the = + projection of the intrinsic state. 

Also in 16 0, the angular motion of an a cluster at the surface is blocked by other three as. The Pauli blocking 
effect between a clusters causes the angular correlation of the tetrahedral configurations in a compact 4a state. For 
transition from the angular correlated cluster state of the ground state to the uncorrelated cluster state like a cluster 
gas, at least two a clusters need to spatially develop to move freely in the allowed region without feeling the Pauli 
blocking effect. The 4a-cluster gas state in 16 has been suggested near the 4a threshold energy. The suggested 
excitation energy E x ~ 15 MeV is almost twice of E x = 7.66 MeV for the 3a-cluster gas state in 12 C. This might 
correspond to the energy cost of the spatial development for two as in 16 compared with that for an a in 12 C. 



(a) (b) 

(y,z) 




FIG. 5: (a) The Pauli allowed and forbidden regions for the rotation of the angle 9 y for the third a in the (y, z)-plane around 
the 2a core. The norm Ns, a {y = r cos 6 y ,z — rsin6 y ) normalized at (y = 0, z — r) is presented. The Nsa ~ 1 an d Ms a ~ 
regions are the Pauli allowed and forbidden regions, respectively, (b) A schematic figure for the position of the third a around 
the 2a core, (c) A schematic figure of the strong cluster coupling state for a compact 3a state with the strong the Pauli blocking 
effect, and (d) that of the weak cluster coupling state. See the text. 
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IV. CLUSTERS ON A FERMI GAS CORE IN ONE DIMENSIONAL BOX 



A. Concept of one-dimensional cluster model 

As mentioned, the triangle deformation of the 12 C ground state can be interpreted as the density wave on the 
axial symmetric oblate state and it is characterized by the static surface density oscillation with the wave number 
three on the edge of the oblate intrinsic state. The symmetry breaking in the ground state originates in the 4-body 
correlation in a clusters. The axial symmetry is broken by the angular correlation between a clusters due to the Pauli 
blocking effect. As a clusters develop, the Pauli blocking weakens and the angular correlation vanishes. Then the 
system transits from the symmetry broken state of the 12 C(0^) to the symmetry restored state of the 12 C(0j), where 
a clusters are freely moving in a dilute density like a gas. 

In this section, we consider Pauli blocking effects in a schematic model of two clusters on a Fermi gas core in a 
one-dimensional (ID) finite box, and discuss how the translational invariance is broken to form an oscillating density 
(inhomogeneous) state and how the broken symmetry is restored to a uniform density (homogeneous) state. 

Let us consider a schematic model for a simplified 3a system consisting of two a clusters around the a core. To 
concentrate only on angular correlations between two as, we ignore the radial coordinate degree of freedom and 
assume that two a clusters are moving at a certain distance r from the core. Taking the body-fixed frame so as one 
of two a clusters on the z-axis and the other a in the y-z plane, then, the angular motion between two as can be 
reduced to the ID problem, where the first cluster sits around the origin and the second cluster exists in a finite size 
L = 2irro box with the periodic boundary condition. For the core effect, we take into account only the Pauli blocking 
from the core particles, which is treated here by a Fermi gas core for simplicity. 

We extend the ID-cluster model and treat the cluster size b, the box size L, and the core Fermi momentum k c 
as parameters. In this model, we do not perform energy variation nor calculate energy eigenstates. Moreover, the 
mechanism of cluster formation and the dynamical change of the cluster structure or the cluster size are beyond our 
scope. We give a model wave function with fixed parameters 6,L,fc c based on simple ansatz, and analyze Pauli blocking 
effect between a cluster and the core and that between two clusters in the given wave function. From behavior of 
the wave function and features of the Pauli blocking effects, we conjecture how the transition between the Fermi gas, 
BCS-like, DW-like, and BEC-like states occurs. 



B. Formulation of one-dimensional cluster model 



We explain details of the ID-cluster model wave function. Two n-body clusters are formed on a Fermi gas core in 
a ID box with the box size L = 2irr . For an a cluster n = 4 and a cluster consists of p f, p I, n f, and n J,. Basis 
single-particle wave functions in the periodic boundary condition are given as e +ikx /V2L and e~ ikx /V2L, where the 
momentum k = lixkjL. Here k is the dimensionless momentum k = krg and takes k = 0, 1, 2, 3, • • • . Core nucleons 
are assumed to form a Fermi surface at k c . It means that k < k c orbits arc occupied by core nucleons and they are 
forbidden single-particle states for constituent nucleons of clusters. Using the dimensionless parameter k c = k c r$, the 
lowest allowed momentum is fci? = kp/ro for Hf = k c + 1 and k > Uf orbits arc allowed. 

We assume the first cluster (ai) wave function that localized around x = in a Gaussian form, 

ip ai = U x (j) (x lx )X x , (9) 
X = T<7 = {pt,pl,nt,nl}, (10) 



O (x) = E /( fc ) ( e+lkx + e ~ lkx ) ■ ( n ) 

VAL k>k F 



f(k) is the Fourier transformation of the Gaussian with the cluster size b and given as 



where no is determined by the condition, 



/(fc)=noexp(-^), (12) 



E i/( fc )i 2 = x > ( 13 ) 

k>k F 



to normalize 4>o(x) in the box. The wave function indicates that four species of nucleons, p f,p i, n 1 s , n i, occupy the 
same spatial orbit with the Gaussian form forbidden partially by the Fermi gas core, lib/ L and kp are small enough, 
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the wave function is localized well around x = 0. Since there is no identical Fcrmion in a cluster, the antisymmctrizer 
can be omitted in the expression. 

Next we consider the second cluster (0:2)- Assuming that «2 is localized around a position x = s, the wave function 
is given by the shifted function, 

r a2 = HUx2 X )X X , (14) 

a 

tj>.{x) = -j= J2 f(k) (e+^-s) + e - ik ^) . (15) 



2L k>k F 

Then, the normalized wave function of the 2a system with the parameter s is written as 

*.=^ B a W^{^^}. ( 16 ) 
Npb(s) = \-^A{foM\ 2 (17) 

= 2 1^0(^1)^8(^2) - Mx2)4's(xi)\ 2 (18) 

= l-((/»o0 s |0^o). (19) 

A/"pb(s) is the overlap norm for two identical nucleons. It is a function of the parameter s for the localization center 
of the second cluster. Afp B means the overlap norm for the 2a system and it is an indicator to evaluate the Pauli 
blocking effect between two clusters. In the case that two clusters feel no Pauli blocking, A/"p B (s) = 1, while in the 
case that they feel complete Pauli blocking, Afp B = 0. It means that A/"p B (s) stands for the Pauli " allowedness" for 
ai2 around s. 

The density distribution for X = tcf particles in the 2a state is 

pX(x) = {* s \^5(x-x i )\V s ). (20) 

Note that the density distributions for all kinds of nucleons are the same in the present cluster model and the total 
nuclear density is ^ x p*{x) = 4p*(x). In a similar way, the x = rcr density for the one-a state ip ai is 

/&(*)H^l£*(*-*0hki>- ( 21 ) 

In the present model, we can express all the formulation with the dimensionless parameters L = L/ro = 2ir, 
b = b/rg, s = s/tq, x = x/ro, k = kro, k c = k c ro, kp = kpr®, and so on. The dimensionless densities are also defined 
as p = pro. Then, the state l^s) is specified with the dimensionless parameters (kF,b, s) and does not depend on 
the scaling factor r^. Since the number of clusters is fixed to be two, a larger volume size L corresponds to a lower 
cluster density. Because of the scaling, a larger volume size L = 2-ktq with a fixed cluster size b is equivalent to a 
smaller cluster size b with a fixed L. That means, the parameter b indicates the cluster size relative to the box size 
and also corresponds to the cluster density. The kp = k c + 1 is the lowest allowed momentum just above the core 
Fermi momentum k c , and it relates to the density of the core particles. 



C. Pauli blocking effect from the core to one cluster 

Let us describe the structure change of one cluster on the Fermi gas core in the ID-cluster model. The coefficients 
f(k) of the Fourier transformation and the density distribution p* (i) for the a\ cluster are shown in Figs. [6] and [7] for 
kp = 1 and kp = 2, respectively. Because of the Pauli blocking effect from the core as well as the finite volume effect, 
the structure of the cluster changes from the original Gaussian form depending of the parameters b and kp. In the case 
of small b, which corresponds to the small cluster size b or the large volume size L, the coefficient f(k) is distributed 
widely and has a long tail toward the high momentum region, and the density is localized well around x = (and 
also x = 2-7T for the periodic boundary). With increase of kp, the low momentum components are truncated and the 
density shows a oscillating tail. Nevertheless, if the cluster size b is small enough, the density is still localized. With 
increase of the cluster size b, the density localization weakens and the density approaches the periodic one. Then, 
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the component of the lowest allowed orbit kp = k c + 1 is dominant and the components of higher orbits decrease. 
The localization declines more rapidly in the case of the larger kp. It means that, as the b increases or as the kp 
increases, the spatial correlation between nuclcons (inner correlation) in a cluster becomes weak. We call the case of 
the weak localization that the wave function has the dominant fcp component and minor k > kp components "the 
weak coupling regime", and the opposite case of the strong localization with significant k > kp components "the 
strong coupling regime". In the weak coupling limit, the one-cluster density goes to cos 2 (fc_Fi)/7r. 



For two a clusters expressed by the cluster wave function \l/ s , the clusters U\ and <xi are assumed to be localized 
around x = and x = s, respectively. The Pauli blocking effect between ot\ and «2 is evaluated by the overlap norm 
A/"p B (.s), which is an indicator for Pauli allowcdness. The Pauli allowedncss Afp B (s) ~ means that the s region is 
blocked by the a± cluster and is a forbidden area for a-i- The middle panels of Figs. [5] and [7] show the allowedncss 
A/p B (s) (thin solid lines) as well as A<pb(s) (dashed lines) plotted as a function of s. In principle, the Pauli blocking 
effect reflects the 1 probability of the ai cluster, and therefore, the forbidden region corresponds to the relatively high 
density region of the a\ cluster. In case of a small cluster size b, the allowed region exists widely and the forbidden 
region exists only in the small area close to x = and i = 2ir. That is to say, the second cluster feels almost no 
Pauli blocking effect except for the region near the first cluster (ai). With increase of b, the forbidden area spreads 
in a wide region around x = and i = 2ix, and the allowed region for ct2 becomes narrow. In the weak coupling 
regime, where the ot\ density is periodic, the allowed region with -<Vp B (s) ~ 1 also shows the periodicity. Reflecting 
the cos 2 (fc^x) periodicity of p ai (x), the areas s w ir2m/2kp (m = 0, ■ • • , 2kp) are the forbidden regions, while the 
areas s = Sj = n(2j — l)/2kp (j = 1, • • ■ , 2kp) are the allowed regions. 

As mentioned, the parameter b corresponds to the cluster size relative to the box size. When the cluster density is 
low enough, the Pauli blocking effect is weak and almost all region is allowed for the ct2 cluster except for the position 
close to the a±. On the other hand, in the case of the high cluster density the Pauli blocking effect is strong and the 
allowed s region is restricted in the periodic regions. 



We first discuss the features of the two cluster wave function ty s with a fixed parameter s. It means that the center 
of the second cluster is located around the fixed position. Later, we will discuss how the spatial correlations between 
clusters (inter-cluster correlation) can be affected by the Pauli blocking. 

We choose Sj — n(2j — l)/2/cp with j = 1 and j = kp, which are the allowed positions Sj nearest to x = 
and x = 7T and corresponds to the smallest and largest inter-cluster (a\-a2) distances, respectively. The density 
distribution p%(x) in the 2a-cluster wave function ^ s are shown in the right panels of Figs. [5] and [7J for kp — 1 and 
kp = 2. In the strong coupling regime, for instance, the (kp, b) = (1.0, 0.25) state, the density shows the clear two peak 
structure and indicates that two clusters are well isolated without almost no overlap. As the b increases, the overlap 
region between clusters gradually increases and the density changes to the oscillation structure, in particular, in the 
case of j = kp. The density oscillation is remarkable, for instance, in the (hp, b) = (1.0, 1.0) and (kp, b) — (2.0, 0.75) 
states, which shows the 2kp + 1 periodicity. With further increase of b, the density oscillation weakens and finally 
disappears to the uniform density, and the system goes to the Fermi gas limit with the Fermi momentum kp. 

In the present model, we put the «2 cluster around the fixed position s. For more realistic wave functions of two a 
clusters, one should extend the model by taking into account the motion of the ai cluster relative to the o,\ cluster. 
Microscopically, it corresponds to superposing ijj^ with various values of the parameter s as 



The spatial correlation between clusters (inter-cluster correlation) is expressed by the weight function F(s). On the 
other hand, the correlation between nucleons inside a cluster (inner correlation) is described by the intrinsic structure 
of a single-cluster wave function of tp ai or ^ s determined by the parameters b and kp, and it is given by hand in the 
present model, where the cluster formation and its intrinsic structure are a priori assumed. 

Moreover, the wave function should be projected to the total momentum Kq = for the center of mass motion 
(c.m.m.) of two clusters so that the translational invariance is restored in the finite system. 



D. Pauli blocking between two clusters 



E. Transitions from strong coupling to weak coupling regimes 




(22) 
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In the following discussions, we do not treat the ci2-cluster motion explicitly, but conjecture how the Pauli blocking 
may restrict the motion of the «2 cluster and affect to the spatial correlations between clusters (inter-cluster corre- 
lations). Since the area of high a\ density is blocked as mentioned, the ai cluster may move in the allowed regions 
with large A/p B (s). We assume that the interaction between clusters is weak and the Pauli blocking effect, which acts 
like an effective repulsion, gives the most important contribution in the relative motion between clusters. 

1. BEC-like state 

Let us consider the strong coupling regime, where the cluster size b is small enough. The A/p B (s) curve shows 
a wide open window of the allowed region as in the case of (kp,b) = (1.0,0.25). Since the Pauli blocking effect is 
weak, the a 2 can move in the wide allowed region almost freely. It means the strong inner correlation but almost no 
inter-cluster correlation. In such the uncorrelated cluster state, two clusters may condensate approximately the zero 
momentum state in the ground state similarly to the BEC phase. 

2. DW-hke state 

As the b increases, the open window for the allowed region closes and a 2 no longer can move freely. Instead, 
the allowed region becomes discrete, and the a2-cluster center may be confined around the allowed s values, Sj = 
7r(2j-l)/2fc F (j = !,■■■ ,2k F ). 

For a possible wave function for 012, one may consider a superposition of wave functions with Sj as 

^ 2 = E^)^i- ( 23 ) 

3 

As seen in the density distribution p%(x) of two cluster states shown in Figs. [5] and [7J the density oscillation shows 
the 2kp + 1 periodicity for any Sj values. 

The density oscillation can be remarkable provided that the amplitude F(sj) for j = kp (nearest Sj to the middle 
point x = 7r) is relatively larger than those for Sj around s — and 2tt because of the effective repulsion between 
clusters due to the Pauli blocking effect. In other words, possible localization of the amplitude F(5j) because of the 
Pauli blocking effect causes the spatial correlation between clusters (inter-cluster correlation). In such the case, the 
density oscillation shows the clear 2kp + 1 periodicity whose origin is the DW-type correlations. Indeed, the state 
contains dominantly the coherent lp-lh components of a ±kp + 1 particle and a ±kp hole on the Fermi surface at 
kp. The lp-lh correlation carries the 2kp + 1 momentum and causes the 2kp + 1 periodicity. This correlation is the 
similar to that of the DW phase. In the weak coupling approximation, the spatial wave function Eq. (|11[) for an a 
cluster is approximated by the dominant fcp component with a minor mixing of the kp + 1 component as 

4>o(x) = (cos kpx + ecos(kp + l)x). (24) 

V 2L 

In this approximation, it can be proved that the 2a wave function *5? s for j = kp actually contains the dominant 
DW-type lp-lh correlation in the particle-hole representation (See appendix IB"]) . 

In the opposite case that there exists attractive force between clusters, the amplitudes F(Sj) may gather to smaller 
Sj. It corresponds to the exciton(Exc)-type correlations described by the coherent lp-lh components of a ±kp + 1 
particle and a ^fkp hole on the Fermi surface as described in appendix[Bj In this case, the lp-lh carries the momentum 
\kp + 1 =p kp\ = 1, which suggests that spatial density oscillation is not so remarkable. Indeed, such the feature is 
seen in the weaker density oscillation in p*(x) for s = Sj (j = 1) shown in Fig. [JJ (dashed lines in the right panels) for 
kp = 2. We should comment that the system for kp — 1 is a special case that the small s region, for instance, the 
region s < 7r/4 is forbidden and the Exc-type correlations are suppressed. 

3. BCS-like state 

With further increase of 6, the kp component becomes more dominant. In the case f{k) <C /(fcf) for k > kp, the 
structure of the Pauli allowed area for a 2 approaches the pure periodic one following the periodicity of the a% density 
Pai(x) ~ cos 2 (kpx)/y/w. Then, a linear combination of ipa 2 with an equal weight may be the lowest state to restore 
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the translational invariance because $ s for different may degenerate energetically. It means that the state no longer 
has the spatial correlation between clusters (inter-cluster correlation) and it corresponds to a BCS-like state. Namely, 
the BCS-like state has the weak inner correlation and no inter-cluster correlation. 

As described in appendix [B] in the weak coupling approximation, the total c.m.m. momentum Kq = state 
projected from the equal weight linear combination of two cluster wave functions s is equivalent to the BCS-like 
state containing 2p-2h configurations of a (/cf + 1, —kp — 1) XaXP particle pair and a (fcp, — kp) XaXP h°l e pair [see 
Eq. (|C5[) ]. In the 2p-2h state, all kinds of pairing is coherently mixed so as to keep the spin-isospin symmetry of a 
clusters. 

4- Fermi-gas state and correlations 

In the large b limit, excited components of k > kp vanish and the system goes to the Fermi gas (FG) state with 
the Fermi surface at kp. Needless to say, the FG state has no inner correlation nor inter-cluster correlation. This is 
nothing but the uncorrelated state. On the other hand, the correlated states are characterized by configurations of 
excited k > kp components. In the weak coupling regime, we can recognize DW-like, Exc-like, BCS-like states, or 
mixing of them by correlations in lp-lh and 2p-2h configurations on the Fermi surface. 

5. Diagram of structure transitions 

As discussed above, in the present schematic model for two clusters on the Fermi gas core in the ID box, the cluster 
state is expected to show the BEC-like, the DW/Exc-likc, BCS-like, or FG behaviors depending on the cluster size b 
and the lowest allowed momentum kp. We here assume the criterion to judge a correlation type for a wave function 
with given b and kp values, and show a diagram of structure transitions on the b-kp plane. 

For the criterion, we define Afc by the deviation of k from the lowest allowed momentum kp. 

( A ^) 2 = E fim-kpf. (25) 

k>k F 

It indicates the deviation from the FG state \HF). In the case Ak <C 1, k > ftp + 2 components are negligible and 
the coefficient e = f(kp + 1) of the kp + 1 component approximates e « Afc. In the Ak = limit, the system goes 
to the FG state. For the criterion to classify a wave function with given b and kp to a correlation type, we use the 
deviation Afc and the Pauli allowedncss N PB (s) at the middle point of the box s = L/2 = tt. 

In Table HI we list the criterion for various correlation types. In the table, we also show the typical examples of the b 
values for kp = 1 and kp = 2 that satisfy the criterion. The densities and the Pauli allowedncss for the corresponding 
states are already shown in Figs. [5] and [7J The BEC-like state is expected to appear when the Pauli allowed region is 
widely open. Then, we adopt the condition J\fp B (s = tt) > 0.8 for the BEC-like state. The DW-like and/or Exc-like 
states may appear when the Pauli allowed region is restricted. For this condition, we use Afp B (s = tt) < 0.1. The 
DW/Exc-likc states may change to the BCS-like state when the ai-cluster density p ai (x) approaches cos 2 (kpx)/ \[tt 
and the Pauli allowed area for a2 becomes almost periodic. For the criterion that p ai (x) cos 2 {kpx) / y/rr , we adopt 
the ratio of the a\ density at x = to that at x = L/2 = tt. In the weak coupling regime of Afc <C 1, the ratio 
Pa ± ( x = 0)/Pqi(^ = 7r ) i s approximately given by 1 — 4e m 1 — 4Afc, therefore, we use Afc for the measure. We apply 
Afc > 0.05 for the DW/Exc-like states. This approximately corresponds to the ratio p&^x = 0)/p* 1 (x = tt) < 0.8. 
With the decrease of Afc, the 2a wave function may gradually change to the FG state via the BCS-like state. Since 
the higher momentum components k > kp nearly equals to (Afc) 2 in the weak coupling regime, we use Afc as the 
measure for structure transitions from the DW/Exc-like to the FG state as listed in the table. For instance, the 
condition Afc < 0.001 for the FG state indicates the contamination of fc > kp components is less than 10~ 6 . 

The diagram of the structure transitions between the FG state, BCS-like state, DW/Exc-like state, and BEC-like 
state on the kp-b plane is shown in Fig. [8j 

For a system with higher kp, nucleons in a cluster couple more weakly to each other because of Pauli blocking 
from core nucleons, and therefore the FG region is wider in the diagram. However, one should care about that the 
assumption of the sharp surface for the core Fermi momentum might be inadequate, in particular, in case of high 
kp. The core Fermi surface may diffuses in correlated states. If the surface diffuses, the lower orbitals below kp arc 
partially allowed for nucleons in clusters. Then, the weakening of the cluster localization by the Pauli blocking from 
core nucleons can be quenched. The present model should be extended by incorporating the surface diffuscness of the 
core Fermi surface. 
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TABLE I: For the criterion to classify the ID-cluster wave function with given 6 and k F into various types of correlation. The 
deviation Afc and the Pauli allowedness Mps{s = L/2 = tt) are used for the criterion. The typical examples of the 6 values for 
fcp = 1 and k F = 2 that satisfy the criterion are also shown in the table. 



correlation criterion example 

k F = 1 k F = 2 

FG Afc < 0.001 6 = 2.0 
FG-BCS crossover 0.001 < Afc < 0.005 6 = 2.0 6 = 1.5 
BCS-like 0.005 < Afc < 0.01 
BCS-DW/Exc crossover 0.01 < Afc < 0.05 6 = 1.5 6 = 1.25 
DW/Exc-like 0.05 < Afc and N PB {L/2) < 0.1 6=1.0 6 = 0.75 
DW/Exc-BEC crossover 0.1 < N PB {L/2) < 0.8 6 = 0.5 6 = 0.5 
BEC-like 0.8 <M *% B {L/2) 6 = 0.25 6 = 0.25 



We also should comment that, in a real system, two parameters b and kp are uncontrollable in general. As explained 
before, b indicates the cluster size relative to the box size (the cluster density), while kF relates to the core density. 
The size of clusters on the Fermi gas core should be determined dynamically as a consequence of nuclear interactions 
between constituent nuclcons in clusters. It should be also affected by the neighboring cluster, i.e., cluster density 
as well as kp. Moreover, the ratio of cluster nucleons to the core nucleons should be determined dynamically. The 
extension of the present model by taking into account dynamical change of cluster structure with the use of effective 
nuclear interactions is an important issue to be solved in future study. 



F. Correspondence to finite nuclei 

As shown in the AMD+VAP calculation, three a clusters are formed in the ground state of 12 C even though 
the existence of any clusters is not a priori assumed in the framework. Once three a clusters are formed, it can 
be associated with the schematic ID-cluster model. Angular correlation between two a clusters around the a core 
corresponds to the spatial correlation between two a clusters in the ID-cluster model with kF = 1- The cluster size 
b is 1.62 fm for v = 0.19 fm~ 2 used in the AMD calculations. If we adopt the r.m.s. matter radius r a = 1.72 fm 
of the core a as the radial size ro, the dimensionless b = b/ro is estimated to be b = 0.94. Or if we use the r.m.s. 
radius of cluster positions of 3 as evaluated from r 2 + 7- 2 = i? 2 , we get b = 0.87. In both cases, b ~ 1. As already 
described, the state with = (1.0,1.0) in the ID-cluster model shows the remarkable density oscillation with 

the wave number three in the DW-like regime. It is consistent with the intrinsic structure with the 12 C(0^ ) obtained 
with the AMD+VAP calculation. Indeed, the 12 C(0+) has the (V 3 ~ 3 — F 3 +3 )/a/2 component and the surface density 
shown in Fig. [5]is similar to the oscillating density in the ID-cluster state for {kF , b) = (1.0, 1.0) in Fig. [51 

As discussed in Ref. [37[, the a correlation in the pentagon ground state of 28 Si can be interpreted as the density 
wave on the sd-shell oblate state. The 28 Si ground state is associated with the ID-cluster model with kp = 2 
considering a core consisting of the spherical 16 and four nuclcons in oblate orbits. Then, the pentagon shape can 
be understood again by the DW-like state in the ID-cluster model wave function with the 2fcp + 1 = 5 periodicity. 

In case of the a correlation in the 16 ground state, the tetrahedron shape can not be connected directly to the ID 
problem. However, when we focus only on the {Y^~ 3 — Y^~ 3 )/y/2 component of the intrinsic density in the 16 ground 
state, the density oscillation is characterized by the wave number three periodicity similar to that of the triangle 
shape in 12 C(0^~) and the deformation feature is associated with the DW-type correlation in the ID cluster model as 
in 12 C. 



V. SUMMARY AND OUTLOOK 



We investigated a-cluster correlations in the ground states of 12 C and 16 while focusing on the surface density 
oscillation in the intrinsic states. The intrinsic states of 12 C and 16 obtained by the AMD+VAP method show 
triangle and tetrahedral shapes, respectively, because of the a correlations. The formation of a clusters in these 
states was confirmed in the AMD framework, in which existence of any clusters are not a priori assumed. The 
intrinsic deformations are regarded as spontaneous symmetry breaking of rotational invariance. It was shown that 
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the oscillating surface density in the triangle and tetrahedral shapes is associated with that in DW states caused by 
the instability of Fermi surface with respect to a kind of lp-lh correlations. 

To discuss the symmetry breaking between uniform density states and the oscillating density states, a schematic 
model of a few clusters on a Fermi gas core in a one-dimensional finite box was introduced. In the model analysis, 
we conjecture structure transitions from a Fermi gas state to a DW-like state via a BCS-like state, and to a BEC-likc 
state depending on the cluster size relative to the box size. 

In both analyses of the BB-clustcr model and the schematic ID-cluster model, Pauli blocking effects arc found 
to play an important role in the DW-like state. The breaking of the translational invariance in the DW-like state 
originates in the Pauli blocking effect between clusters, which acts as an effective inter-cluster repulsion and restricts 
cluster motion. 

In the present analysis with the schematic ID-cluster model, we do not perform energy variation nor calculate 
energy eigenstates. Moreover, the mechanism of cluster formation and the dynamical change of the cluster structure 
or the cluster size are beyond our scope in the present paper. The extension of the present model by taking into 
account dynamical change of cluster structure with the use of effective nuclear interactions is an important issue to 
be solved in future study. Also the cluster and core formations as well as diffuseness of the core Fermi surface should 
be studied in more realistic models. Furthermore, the assumption that the interaction between clusters is weak and 
the Pauli blocking effect gives the major contribution to the inter-cluster motion may be too simple. To clarify which 
state of BCS-like, DW-like, or EXc-like ones realizes it is essential to explicitly solve the problem of inter-cluster 
motion by taking into account nuclear forces or inter-cluster interactions. 

It would be interesting to associate the present picture for clusters in the ID finite box with phase transitions in 
infinite matter. In the extension of the present model to infinite matter problem, one should take care of the differences 
between finite systems and infinite systems as follows. Firstly, momentum k values in a finite box is discrete because of 
boundary condition, while they are continuum values in infinite matter. In the description with discrete momentum, 
long range correlations beyond the box size L is not taken into account. What we call the BCS-type correlation in 
the present ID-cluster model is the correlation in the range of the box size at most. In second, the total momentum 
of c.m.m. should be projected to zero in finite system, while it is not necessarily zero in infinite systems. In spite of 
those differences, the ID-cluster model may give a hint to understand an origin of DW in infinite matter. 
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Appendix A: Triangle and tetrahedral deformations in the 3a- and 4a-cluster structures 

We here describe triangle and tetrahedral deformations and particle-hole representations of 3a- and 4a-cluster 
wave functions. As explained in Ref. J37|, the Brink-Bloch (BB) 3a-cluster wave functions [53| with a regular triangle 
configuration in the small inter-cluster distance case can be rewritten as, 

*ta to « II {^oo^(0i-i + £02+2)^(01+1 - e0 2 _ 2 )# TCT } . (Al) 

rcr 

Here </>; m is the single-particle orbit in the harmonic oscillator potential, X Ta is the spin-isospin wave function with 
t = {p,n} and a = {tii}- e is the order parameter for the axial symmetry breaking and it corresponds to u/v 
in the BCS theory. Here e is considered to be small to omit e 2 and higher terms and taken to be a real value. In 
principle, e can be a complex value but the phase is arbitrary and can be changed by the rotation around the z axis 
(transformation of the polar coordinate <p). The density of the $£^(e) state is 

= "W e-fl \ l + 2 i Sin2 ^ + e2 ^S sin3 (0)( e " 3 ' V - + °(t 2 )) > (A2) 
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FIG. 6: (Color on-line) The results of one and two clusters on the Fermi gas core in a one-dimension box. Left: The coefficients 
f(k) of the Fourier transformation. Middle: The density Pq 1 (5) of the first a located around x = (red thin lines), the PB 
effect Afps(s) for the second a (black solid lines). ATpb{§) is also shown (blue dashed lines). Right: The density p* (x) (black 
solid lines). The s for the center position of ol-z is chosen to be §j with j — kF, which is the closest §j to x = L/2 = 7r among 
the allowed s values. The density p^ 1 (i) is also shown for comparison(red thin lines). The results for the cluster size 6 = 2.0, 
1.5, 1.25, 1.0, 0.75, 0.5, 0.25 are shown. 



and its multipolc decomposition at a certain radius r = Rq is 
p{R ,6,ip) 



r 3 (Y-\e^) Y+ 3 (e,<p) 



0(e 2 



(A3) 



V2 V2 

In a similar way, the BB 4a-cluster wave functions with a regular tetrahedral configuration can be rewritten as 

*£f (e) ~ II WooAW^ioAV^i-i + e02+i)AV CT (0i+i + e<t> 2 -i)X T <r} ■ (A4) 



Here it is assumed that the inter-cluster distance in the tetrahedral configuration is small, i.e., e is small. The density 
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FIG. 7: Same as Fig. [6] but for Hf = 2. The blue dashed lines indicate the the density p* (x) for the smallest allowed s value, 
Sj with j — 1. The results for the cluster size b — 2.0, 1.5, 1.25, 1.0, 0.75, 0.5, 0.25 are shown. 



P^) = ^m& e ~^ { 1 + 2 ^ sin2 ^ + sin 2 W(e 2 ^ + e~ 2 ^) + 0( £ 2 )| , (A5) 



and its multipole decomposition at a radius r = Rq is 



= (i + ifVWM + «V iiifl-^Tf + vf^ ■ (A6) 



r 3 " 2 (6», ip)/V2 + Y+ 2 (6, <p) /y/2 in the e term can be also transformed to 

^-Y*{6, <p) + ^-Yi\e, <p) - ^-Y+ 6 {G, ip), (A7) 



^2 V2 3 3 v ,yj 3 3 v ^' 3 3 
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by a f2 rotation Rq. 

By using the creation and annihilation operators, aj m and ai m . x for the <\>i m X T(J state with x = TtT ; the $£^(e) 
and &fa( e ) states are expressed as 

l^fcf(e)) = n a oo,x( ai - 1 .x + e 4+2, x )(ai+i, x -e4-2, x )|-) ! (A8) 
x 

Ma (e)> = II a ™,x at io,x( a i-i.x + a4+i, x )(oi+i,x + «4-i, x )|->- (A9) 

X 

In the particle and hole representation, they are rewritten as 

X 

\0) F = ]J (oSo )X ai_i lX oi+i lX ) h), (All) 

X 

and 

l*£f W) = + e 4 + i,x 6t i + i,x)( 1 + e 4-i, x &i-i,x)l°>^ ( A12 ) 

X 

|0)f = (4o,x a io,x°i-i.x a i+i,x) I")- ( Al3 ) 



Here the hole operator bj m is defined to be b\ m = a;_,„ iX . Equations (|A10|) and (|A12[) indicate that the 3a- and 
4a-cluster wave functions contains the DW-type lp-lh correlations carrying finite angular momenta. 



Appendix B: Weak coupling regime of one-dimensional cluster model 

We consider the weak-coupling case of the ID-cluster model described in llVI Hereafter, we take ro = 1 and consider 
the case that all the parameters equal to the dimensionless parameters, for instance, b = b, k = k, and s = s. We 
define the creation and annihilation operators a k and a^ fcx for the single-particle states <pkX x and <j>kX x , where 

<f)k and 4>-k are momentum fc and — k plane waves written as 4>k(%) = e +lkx /\/2L and <f>-k{x) = e~ lkx /\/2L in the 
coordinate representation. The core state \C) where all k < kp states are occupied by core nucleons is written as 

1^=11 n «la«i->. (bi) 

x |fc|<fc F 

where |— ) is the vacuum. 
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In the case of weak coupling, the coefficient f(k) of the momentum k and — k plane waves in Eq. (|9|) is dominant 
for k = k F and it rapidly decreases as k increases. In the limit, \f(kp)\ ^> \frkp + 1)| >• \frkp + 2)| • • • , the leading 
term of the 2a-cluster state in the ID-cluster model is the Fermi gas state with the surface momentum kp, which is 
the uncorrelated state. For the correlated system, the deviation from the Fermi gas state is important. We take into 
account the next leading term, k = kp + 1 components, and ignore the higher momentum components for k > kp + 2. 
We set f(kp) = 1 and f(kp + 1) = e with e << 1 and take the leading term of e. Then the state consisting of two a 
clusters around x — and x — s written by Eqs. (jH]) and (|14[) on the Fermi gas core is expressed as 

\^ ai r a2 C)^l[{(c*~co)al x al kFX 

x 

+ e(c - Cl )a^ x 4 FX + e(c* - c^al^al^ 
+ e{c* a - Cl )al Fx al kFX + e(c - c^al^a^ 

+ ei(c* 1 -c 1 )ai Fx al k , FX + ---}\C), (B2) 



where k' F = kp + 1, cq(s) = e~ lkpS , and ci(s) = e~ lkpS . The first term is the Fermi gas state with the Fermi surface 
at kp. Its coefficient Cq — cq = 2isin(kps) vanishes for s = irm/kp (m is an integer) because of Pauli principle. In 
other words, the area around s = it2m/2kp (m = 0, • • • , 2kp) is forbidden for the oti cluster center position, while 
the area around Sj = ir(2j — l)/2kp (J = 1, ■ ■ ■ , 2k p) is allowed. 

In the particle- hole representation with respect to the Fermi surface at kp, the second and third terms correspond 
to the Ip — lh states with the DW-type correlation with the wave number k' F + kp = 2kp + 1, and the forth and 
fifth terms are those with the Exc-type correlation with the wave number k' F — kp = 1, and the last term is the 2p-2h 
state. We rewrite the coefficients in \ip ai ipa 3 C) as 



k F x -k F X 
X 



+ cC D w{s)al FX a{ FX + eC* DW (s)a f _ k , FX a f _ kFX 
+ eC Ex c(s)al FX al kpx + eC* Exc {s)a}_ k , FX al FX 

+ ^C 2 (s)al FX al k , FX + ---}\C), (B3) 

(B4) 

(B5) 



C (s) = 


Co 


- c 




— 2i s'm(kps), 


Cdw(s) 




co - 


Cl 


e~ ls - 1 




r* — 
c 


Co 


2 ' 


C Exc {s) 




c o - 


Cl 


e~ ls + 1 




r* — 


CO 


2 ' 


C 2 (s)^ 


4 


- Cl 

- c 




cos(js). 



(B6) 
(B7) 



By using the hole operator b\ = a-fe x we obtain 

\^ i r a2 C)=nlY[C Q (s){l 



\ 



- t C Dw(s)al, FX bi FX + eC* DW (s)al k , FX bl kFX 
+ eC Exc ( S )a{, FX bl kFX - eC* Exc (s)al k , FX bt FX 

+ ^C2( S )al FX al k , FX bl kFX b{ FX + ■■■} \HF), (B8) 

where the Hartree-Fock vacuum \HF) = Il x a I FX a -fe FX |C)- When s is the allowed Sj close to L/2 = ir, i.e., the 
ct2 center is located around the middle of the box, e~ ls approaches —1, and therefore, the DW correlation becomes 
dominant and the Exc correlation is minor, while they are opposite for Sj close to 0. In both cases, the correlation 
disappears in the e — > limit, and the system goes to the Fermi gas limit. 
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Appendix C: Transition between DW/Exc-like and BCS-like states 

We extend the model in the weak coupling limit described in the previous section by superposing all allowed Sj 
states for the ct2 wave function, 

i*} = E F ( s ^i^ c >- ( C2 ) 

3 

In principle, one should determine the coefficients F(sj) following the energy variation to properly take into account 
dynamics of relative motion between clusters, a% and ai. However, we here do not perform energy optimization but 
adopt a simple ansatz for F(sj) as follows. 

As described before, if the Sj ~ L/2 (sj ~ 0) components are main, the DW (Exc) correlation is dominant. In more 
general, in case that the effective interaction between a clusters is repulsive, the amplitudes may be relatively smaller 
in the region around Sj ~ (and also around Sj ~ L due to the periodic boundary) than that around Sj = L/2 = it 
resulting in the DW dominance. On the other hand, if the interaction between a clusters is attractive, the amplitudes 
concentrate on the region around Sj ~ 0, which may cause the Exc dominance. In both cases, spatial correlation 
between ct\ and ui remains and the system shows density oscillation and the symmetry breaking of translational 
invariance. 

Next we construct the state where there is no spatial correlation between a,\ and a,i (no inter-cluster correlation) 
and the translational invariance is completely restored by the superposition and the Kq = (total momentum of 
c.m.m.) projection of the \ip ai ipa 2 C} . The state with no spacial correlation between cl\ and a-i can be constructed 
by the projection to the relative momentum k ai — k a . 2 — state. Here k ai 2 stands for the momentum of c.m.m. of 
011,2' In the weak coupling regime, this is performed by superposing all allowed Sj states with an equal weight as 



3 

Taking averages of Cdw(J)i CexcU) ano - ^(j) with respect to j, we obtain 

\M a2 C) cx I] {l - |(at, Ax - °Ux b -*«) 

x 

--(at b\ -a f ,, bt ) + ■■■ )\HF), (C4) 

This corresponds to the intrinsic state for no correlation between clusters. In the ground state of the finite system, the 
total momentum Kq = k ai + k a2 should be zero. With the projection onto Kq = 0, wc finally obtain the symmetry 
restored state, 

P KG=0 \^a,ipa 2 C) oc J 1 + — V (at a f ,, -at a 1 ,, )(bl b\ -b\ b\ ) + ■■■ \ \HF). (C5) 

\ XcXfi ) 

It is found that the state contains 2p-2h configurations and it is consistent with a BCS-like state, where two particles 
form a (k' F ,— k' F ) XaXp P arr an d two holes do a (Hf,— kp) XaXp pair. In the 2p-2h state, all kinds of \ pairing is 
coherently mixed so as to keep the spin-isospin symmetry originating in the symmetry of a clusters. 

Note that, in the inter-cluster correlated case of the DW dominance state, the Kq = projected state contains 
a] , a b\ b , + at a u bt b , in the e 2 term, which means that lp-lh of y„ always carries 

k'pXc, -k' F Xfi k F x a -k F XH k' F Xf> -k' F Xc k FXf) -kpXc ' 

the finite momentum ±(k' F + &f) as a result of spatial correlation between clusters. 
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